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We recapitulate on aspects of Dyson- Schwinger equation studies relevant to pseu- 
i doscalar mesons: lattice confirmation of the DSE prediction that propagators are 

I nonperturbatively dressed in the infrared; and exact results, e.g., in the chiral limit 

the leptonic decay constant vanishes for every pseudoscalar meson except the pion. 
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1. Gap equation 



^ ' The vast body of pion data available provides compelling evidence that this com- 

. posite particle is a Goldstone mode of the strong interaction associated with the 

dynamical breaking of chiral symmetry. Therefore a legitimate understanding of 
pion observables, including its mass, decay constants and form factors, requires 
' that an approach possess a well-defined and valid chiral limit. This is impossible 

. without a detailed grasp of the connection between current- and constituent-quarks. 

' In QCD the running quark mass is obtained from the solution of 

5-i(p) = Z2(z7-p + mbarc)+ ^1 / 9^ D ^^ip ~ q) S (q) T^^iq; p) . (1) 

pH I This is the Dyson-Schwinger equation (DSE) for the dressed-quark self energy 

or, equivalently, QCD's gap equation, and it is a keystone in understanding dy- 
namical chiral symmetry breaking (DCSB) and the relation between current- and 
i constituent-quarks. On the right hand side of Eq. D^i^^{k) is the dressed-gluon 

\ propagator; r°(g;p) is the dressed-quark-gluon vertex; mbarc is the A-dependent 

current-quark bare mass; and :— J d'^q/ {2'it)'^ represents a translationally- 
invariant regularisation of the integral, with A the regularisation mass-scale. In 
addition, Zi^2{Ck^) are the quark-giuon-vertex and quark wave function renor- 
malisation constants, which depend on A and the renormalisation point, C, as does 
the mass renormalisation constant Zjn{C^^) = Zi(C A"^) / Z2{C K^). The solu- 
tion of Eq. Q has the form 

S-\p) = Z7 -pAip^e) + B{p\e) ^ ^(^2 ^2) • P + , (2) 

where M(C^) = m(C) := mbaro(A)/^m(C^, A^) is the nmning quark mass. 
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Figure 1: Left Panel - Dashed-curve: gap equation's solution in the chiral limit; 
solid curves: solutions for A/(p^) obtained using the current-quark masses in Eq. 

(From Ref. ^.) Data, upper three sets: lattice results for M(p^) in GeV at 
am values in Eq. Q); lower points (boxes): linear extrapolation of these results 
to am — 0. Right Panel - Dashed curve, Z{p'^), and solid curve, M{p'^) calculated 
from the gap equation with m{() = 55 MeV @]. Data, quenched lattice-QCD results 
for M{p'^) and Z{p'^) obtained with am = 0.036 {Z{p'^) is dimensionless.) 



The behaviour of the nonperturbative solution of QCD's gap equation is a long- 
standing prediction of DSE studies |1I2I3) . and typical results are illustrated in 
Fig. n One critical feature is that so long as the kernel of the gap equation has 
sufficient integrated strength at infrared momenta, one obtains a nonzero running 
quark mass even in the chiral limit: 

Z4(C,A)m(C) = 0, A»C. (3) 

This effect is DCSB. It is impossible at any finite order of perturbation theory and 
apparent in Fig. ^ 

The dressed-quark propagator can be calculated in lattice-regularised QCD. Re- 
sults are available in the quenched truncation, and depicted in Fig. ^ arc those of 
Ref. |5] obtained with the current-quark masses {( = 19 GeV) 



a miattico 


0.018 


0.036 


0.072 
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0.030 


0.055 


0.110 



(4) 



The precise agreement with DSE results is not accidental. The essential agreement 
between lattice results and DSE predictions was highlighted in Refs. ^ but Ref. 
0] pursued a different goal. Only recently has reliable information about the gap 
equation's kernel at infrared momenta begun to emerge, in the continuum \7\ and on 
the lattice 0. Reference ^ therefore employed an Ansatz for the infrared behaviour 
of the gap equation's kernel in order to demonstrate that it is possible to correlate 
lattice results for the gluon and quark Schwinger functions via QCD's gap equation. 
This required the gap equation's kernel to exhibit infrared enhancement over and 
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above that observed in the gluon propagator alone, which could be attributed to an 
amplification of the dressed-quark-gluon vertex whose magnitude is consistent with 
that observed in quenched lattice estimates of this three-point function [5] . 

2. Hadrons 

It is evident that reliable knowledge of QCD's two-point functions (the propagators 
for QCD's elementary excitations) is available. Direct comparison with experiment 
requires an equally good understanding of bound states. Progress here has required 
the evolution of an understanding of the intimate connection between symmetries 
and DSE truncation schemes. The best known scheme is the weak coupling ex- 
pansion, which reproduces every diagram in perturbation theory. This scheme is 
valuable in the analysis of large momentum transfer phenomena because QCD is 
asymptotically free. However, it precludes any possibility of obtaining nonpertur- 
bative information, and bound state phenomena are intrinsically nonperturbative. 

The properties of the pion are profoundly connected with DCSB, and chiral 
symmetry and its breaking are expressed through the axial-vector Ward-Takahashi 
identity (fc± = k± P/2, {r-', j — 1, 2, 3} are the Pauli matrices) 

P^rf^ (fc; P) = 5- 1 )z75 y *75 y 5- 1 (fc_ ) - iM {() (fc; P) - (fc; P) zM (C) . 

(5) 

This identity connects the axial- vector vertex: rg^(fc;P), P is the total momen- 
tum, with the dressed quark propagator: S = diag[5'„, S'd], the pseudoscalar vertex: 
Fg(fc;P), and the current-quark mass matrix: A4{C) = diag[TOu(C), md(C)]- The 
propagator satisfies the gap equation but the vertices are determined by inhomoge- 
neous Bethe-Salpeter equations; e.g., 



rL(fc;P) 



= Zo 



^[xUq;P)]srKi:{q,k;P), (6) 
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wherein Xz^il^P) •= '5((;+)F|^j(g; P)iS(g_) and K{q,k;P) is the dressed-quark- 
antiquark scattering kernel. The importance of DCSB entails that any truncation 
useful in understanding low energy phenomena must be nonperturbative and pre- 
serve Eq. |SJ), without fine tuning. This nontrivial constraint cannot be satisfied 
without an intimate connection between K(q, k; P) and the gap equation's kernel. 

One systematic truncation scheme has been identified that explicates this con- 
nection and hence preserves QCD's global symmetries ^Ol- It is a dressed-loop 
expansion of the dressed-quark-gluon vertices that appear in the half-amputated 
dressed-quark-antiquark scattering matrix: S^K . The leading order term is the 
renormalisation-group-improved rainbow-ladder truncation, which underlies a one- 
parameter model of the quark-quark interaction used successfully in ah initio cal- 
culations of vector and flavour nonsinglet pseudoscalar meson properties 

The existence of a nonperturbative, systematic and symmetry preserving trun- 
cation scheme enables exact results to be proved. For example, it is a general 
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feature of QCD that the axial- vector and pseudoscalar vertices exhibit poles when- 
ever = —rn^^ , where m,r„ is the mass of the pion or any of its radial excitations.* 
This can be expressed for the axial- vector vertex via 



+ n,ik;P) + j^f^TiJk;P), (7) 

where: Fr, Gr, Hr and f |^ are regular as P^ ^™^„' ^^^^^(fc; P) ^ O(P^) and 
nonsingular; Tl^^{k; P) is the ^ bound state's Bethe-Salpeter amplitude: 

riJk;P) = T^-f5[iE^Jk;P)+-/-PF^,Xk;P) 

+ j-kk- PG^^{k;P) + <j^, k^P, (fc; P)] , (8) 

which is determined by the homogeneous Bcthc-Salpeter equation 

[TiJk;P)]^^ = r[xUq-^P)]srKl^{q,k;P); (9) 

Jq 

and /jr,, is the pseudoscalar meson's leptonic decay constant 

A„ 6'^^ p, = Z2 tr r iT*757Mxi„ (g; p) , (lo) 

Jq 

where the trace is over colour, flavour and spinor indices. Equation H10() is the 
expression in quantum fleld theory for the pseudowector projection of the meson's 
Bethe-Salpeter wave function onto the origin in configuration space. 
The analogous expression in the case of the pseudoscalar vertex is 

'^'s^^'^^Ip^W, «o = regular terms + ^^^-"^^ Jfc; P) , (11) 

^p^.iOS^' = z^tT r ir^j,xUq;P)- (12) 

Jq 

Equation (|12|l expresses the pseudoscaZar projection of the meson's Bethe-Salpeter 
wave function onto the origin in configuration space. 

Inserting Eqs. iQ, ifTTl) into Eq. ©, and subsequently equating pole terms, one 
obtains the modcl-indcpcndcnt result 

U^ml^ = K(C) + maiC)] P.AC) ■ (13) 



*The hadron spectrum exhibits a sequence of J^'^' = 0~+ mesons, with 7r(140) being the 
lowest mass entry. In quantum mechanical models the other members of this sequence: 7r(1300), 
7r(1800), . . . , are described as radial excitations of the 7r(140). Aspects of this interpretation persist 
in Poincare covariant studies in quantum field theory and hence we retain the nomenclature. 
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In the chiral limit, Eq. ||3Jl, the axial- vector Ward-Takahashi identity becomes 

P^n^ik;P)^S-'{k+)zj,^+zj,^S-'{k_). (14) 

Assume that chiral symmetry is dynamically broken so that the dressed-quark prop- 
agator has a nonzero Dirac-scalar term; viz., i? 7^ in Eq. 0. It then follows 
that there is a massless pseudoscalar bound state, ni-^g = 0, for which 

flE^,,{k;0) = B{e), (15) 

with similar relations for the other scalar functions in Eq. (|HJ|. This is clearly the 
ground state and finiteness of the right hand side in Eq. (fT5|) entails f^^ ^ 0. 
Moreover, for the ground state pion in the chiral limit 

P.o(C)-P^o(C)-^^4tr /'%"(g)^-J^(g-q)0, (16) 

where (gg)" is the vacuum quark condensate. Hence the Gell-Mann-Oakes-Renner 
relation for the ground state pion appears as a corollary of Eq. H13() . Another 
important corollary of Eq. (|13|l . valid for pseudoscalar mesons containing at least 
one heavy-quark, is described in Ref. [T^ . 

It is plain that m-n-n^o ^ ™'T^n=a ^'^'^ ^ radial excitations and hence in the chiral 
limit m-^^_^g > 0. In this limit it is impossible to avoid the fact that the absence of 
a pole contribution on the right hand side of Eq. I|14|) forces 

/°„,o=0; (17) 

viz., in the chiral limit the leptonic decay constant vanishes for every one of the 
pion's radial excitations.^ In general 



3. Quantitative illustration 

The manner in which these results are realised in QCD can be illustrated using 
the one-parameter renormalisation-group-improved ladder model for the quark- 
antiquark scattering kernel introduced in Ref. |14| . and reviewed in Ref. P]. One 
first solves the gap equation, Eq. whose solution is required to complete the 
specification of the Bethe-Salpeter equation, Eq. 0, and then solves this for the 
pion and and its first radial excitation. 

discussion of this result in chiral quark models is presented in Ref. 1131 . We thank M.K. 
Volkov and V.L. Yudichev for bringing this to our attention. 
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The homogeneous Bethe-Salpeter equation is an eigenvalue problem, with the 
bound state masses: = ~"^^„ ' being the eigenvalues. We want only the first two 
eigenvalues and eigenvectors. They can be obtained from the modified equation 

liP')[riJk;P)]^^= t[xiS<l-,P)]srKi:{q,k-P), (19) 

Jq 

with 1{P^) a scalar, which has a solution for every value of P^ and can therefore be 
solved by iteration. To explain this, consider the equation written in the form 

1{P^) \g) = M(P2) \g) , (20) 

where the matrix M{P'^) denotes the full kernel of the Bethe-Salpeter equation. 
One fixes a value of P^ and "guesses" a solution: |5(0)). The kernel, M{P'^), 
has a complete set of real eigenvectors \gi) with eigenvalues Ai, ordered such that 
Ao > Ai > . . ., and therefore 

oo 

|g(0)) = ^a,|g,), (21) 

i=0 

where ai are real constants and the vector is canonically normalised. It is clear that 

oo r DO ^ 

M(PTl5(0))-E^f =Ao~ «o|5o) +E (22) 

and hence for sufficiently large N, 

M{P^f+' |.g(0)) « Ao Af(P^)^ |<?(0)). (23) 

Thus repeated operation of the kernel on the initial "guess" produces the largest 
eigenvalue, Zo(P^) = Aq, and its associated eigenvector to any required accuracy. 

One completes this exercise for a range of values of P^ and thereby obtains 
a trajectory lo{P^) that maps the P^-evolution of the integral equation's largest 
eigenvalue. It is then straightforward to find that P^ for which Iq{P'^) — 1. This 
is the solution of Eq. ^ so that P^ = —ni^^ and the associated eigenvector is the 
ground state pion's Bethe-Salpeter amplitude. 

The procedure can also be applied to determine the first excited state. One fixes 
P^ < — wjg, and finds the largest eigenvalue and associated eigenvector for this new 
mass-scale. That completed, one again "guesses" the Bethe-Salpeter amplitude but 
now projects out the eigenvector associated with the largest eigenvalue at this P^: 

l5(0)) = |5(0)) - \go) (.9ol.9(0)) . (24) 

The iterative procedure then applied as before to \g{0)) yields the second largest 
eigenvalue and its associated eigenvector. One thereby obtains the P^-evolution of 
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the second largest eigenvalue, Zi(P^). The solution of li{P^) = 1 gives the mass of 
the first excited state and at this the eigenvector is that state's Bethe-Salpeter 
amplitude. Any finite number of excited states can be studied in this way. 

We have not yet specified the meaning of the inner product used implicitly in 
Eq. ()24|l . For models in the class characterised by the rainbow-ladder truncation 

{h\g) tr T h{q; -P) S{q+) g{q; P) S{q.) . (25) 

Jq 

The condition {h\g) — then expresses a statement that the momentum-dependent 
g h vacuum polarisation (overlap amplitude) vanishes at P^, which is akin to 
saying that the P^-dependent mass-mixing matrix for the states g, h does not 
possess off-diagonal terms. 

We have obtained the mass and amplitude for the ground state pion, using the 
complete expression in Eq. JSJ), and found: /^^ = 0.092 GeV; TOtto = 0.14 GeV; 
Ptto = (0.81 GeV)^, at a current-quark mass r7ic;(lGeV) — r7ii,(lGeV) — 5.5 MeV, 
reproducing the results in Ref. |14| . 

Our study of the first excited state is in its early stages and hitherto we have 
only employed the leading amplitude, -Btti, in Eq. and therewith obtained an 
estimate for the mass: 

m^„^, w l.lGeV. (26) 

Assuming P7r„=i ^ P-woi which is supported by our preliminary estimates, it follows 
from Eq. (|18|l that at the physical current-quark mass 

A„=i < 0.016 = 1.5 MeV. (27) 



4- Epilogue 

We have necessarily been brief. There are many other applications of interest to this 
community, among them the ab initio calculation of electromagnetic and transition 
pion form factors, and a calculation of the pion's valence-quark distribution function 
whose discrepancy with extant data raises difficult questions. These and other 
studies are reviewed in Ref. P] . A pressing contemporary challenge is the extension 
of the framework to the calculation of baryon observables, aspects of which are 
beginning to be understood |15j . 
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